Comprehensive analysis of the microbiome in Apis cerana honey highlights honey as a potential source for the isolation of beneficial bacterial strains.

Background: Honey is a nutritious food made by bees from nectar and sweet deposits of flowering plants and has been used for centuries as a natural remedy for wound healing and other bacterial infections due to its antibacterial properties. Honey contains a diverse community of bacteria, especially probiotic bacteria, that greatly affect the health of bees and their consumers. Therefore, understanding the microorganisms in honey can help to ensure the quality of honey and lead to the identification of potential probiotic bacteria. Methods: Herein, the bacteria community in honey produced by Apis cerana was investigated by applying the next-generation sequencing (NGS) method for the V3-V4 hypervariable regions of the bacterial 16S rRNA gene. In addition, lactic acid bacteria (LAB) in the honey sample were also isolated and screened for in vitro antimicrobial activity. Results: The results showed that the microbiota of A. cerana honey consisted of two major bacterial phyla, Firmicutes (50%; Clostridia, 48.2%) and Proteobacteria (49%; Gammaproteobacteria, 47.7%). Among the 67 identified bacterial genera, the three most predominant genera were beneficial obligate anaerobic bacteria, Lachnospiraceae (48.14%), followed by Gilliamella (26.80%), and Enterobacter (10.16%). Remarkably, among the identified LAB, Lactobacillus kunkeei was found to be the most abundant species. Interestingly, the isolated L. kunkeei strains exhibited antimicrobial activity against some pathogenic bacteria in honeybees, including Klebsiella spp., Escherichia coli, Enterococcus faecalis, Pseudomonas aeruginosa and Staphylococcus aureus. This underscores the potential candidacy of L. kunkeei for developing probiotics for medical use. Taken together, our results provided new insights into the microbiota community in the A. cerana honey in Hanoi, Vietnam, highlighting evidence that honey can be an unexplored source for isolating bacterial strains with potential probiotic applications in honeybees and humans.

Molecular identification of Hymenopteran insects collected by using Malaise traps from Hazarganji Chiltan National Park Quetta, Pakistan.

The order Hymenoptera holds great significance for humans, particularly in tropical and subtropical regions, due to its role as a pollinator of wild and cultivated flowering plants, parasites of destructive insects and honey producers. Despite this importance, limited attention has been given to the genetic diversity and molecular identification of Hymenopteran insects in most protected areas. This study provides insights into the first DNA barcode of Hymenopteran insects collected from Hazarganji Chiltan National Park (HCNP) and contributes to the global reference library of DNA barcodes. A total of 784 insect specimens were collected using Malaise traps, out of which 538 (68.62%) specimens were morphologically identified as Hymenopteran insects. The highest abundance of species of Hymenoptera (133/538, 24.72%) was observed during August and least in November (16/538, 2.97%). Genomic DNA extraction was performed individually from 90/538 (16.73%) morphologically identified specimens using the standard phenol-chloroform method, which were subjected separately to the PCR for their molecular confirmation via the amplification of cytochrome c oxidase subunit 1 (cox1) gene. The BLAST analyses of obtained sequences showed 91.64% to 100% identities with related sequences and clustered phylogenetically with their corresponding sequences that were reported from Australia, Bulgaria, Canada, Finland, Germany, India, Israel, and Pakistan. Additionally, total of 13 barcode index numbers (BINs) were assigned by Barcode of Life Data Systems (BOLD), out of which 12 were un-unique and one was unique (BOLD: AEU1239) which was assigned for Anthidium punctatum. This indicates the potential geographical variation of Hymenopteran population in HCNP. Further comprehensive studies are needed to molecularly confirm the existing insect species in HCNP and evaluate their impacts on the environment, both as beneficial (for example, pollination, honey producers and natural enemies) and detrimental (for example, venomous stings, crop damage, and pathogens transmission).

Comparing modern identification methods for wild bees: Metabarcoding and image-based morphological taxonomic assignment.

With the decline of bee populations worldwide, studies determining current wild bee distributions and diversity are increasingly important. Wild bee identification is often completed by experienced taxonomists or by genetic analysis. The current study was designed to compare two methods of identification including: (1) morphological identification by experienced taxonomists using images of field-collected wild bees and (2) genetic analysis of composite bee legs (multiple taxa) using metabarcoding. Bees were collected from conservation grasslands in eastern Iowa in summer 2019 and identified to the lowest taxonomic unit using both methods. Sanger sequencing of individual wild bee legs was used as a positive control for metabarcoding. Morphological identification of bees using images resulted in 36 unique taxa among 22 genera, and >80% of Bombus specimens were identified to species. Metabarcoding was limited to genus-level assignments among 18 genera but resolved some morphologically similar genera. Metabarcoding did not consistently detect all genera in the composite samples, including kleptoparasitic bees. Sanger sequencing showed similar presence or absence detection results as metabarcoding but provided species-level identifications for cryptic species (i.e., Lasioglossum). Genus-specific detections were more frequent with morphological identification than metabarcoding, but certain genera such as Ceratina and Halictus were identified equally well with metabarcoding and morphology. Genera with proportionately less tissue in a composite sample were less likely to be detected using metabarcoding. Image-based methods were limited by image quality and visible morphological features, while genetic methods were limited by databases, primers, and amplification at target loci. This study shows how an image-based identification method compares with genetic techniques, and how in combination, the methods provide valuable genus- and species-level information for wild bees while preserving tissue for other analyses. These methods could be improved and transferred to a field setting to advance our understanding of wild bee distributions and to expedite conservation research.

Unique spatially and temporary-regulated/sex-specific expression of a long ncRNA, Nb-1, suggesting its pleiotropic functions associated with honey bee lifecycle.

Honey bees are social insects, and each colony member has unique morphological and physiological traits associated with their social tasks. Previously, we identified a long non-coding RNA from honey bees, termed Nb-1, whose expression in the brain decreases associated with the age-polyethism of workers and is detected in some neurosecretory cells and octopaminergic neurons, suggesting its role in the regulation of worker labor transition. Herein, we investigated its spatially and temporary-regulated/sex-specific expression. Nb-1 was expressed as an abundant maternal RNA during oogenesis and embryogenesis in both sexes. In addition, Nb-1 was expressed preferentially in the proliferating neuroblasts of the mushroom bodies (a higher-order center of the insect brain) in the pupal brains, suggesting its role in embryogenesis and mushroom body development. On the contrary, Nb-1 was expressed in a drone-specific manner in the pupal and adult retina, suggesting its role in the drone visual development and/or sense. Subcellular localization of Nb-1 in the brain during development differed depending on the cell type. Considering that Nb-1 is conserved only in Apidae, our findings suggest that Nb-1 potentially has pleiotropic functions in the expression of multiple developmental, behavioral, and physiological traits, which are closely associated with the honey bee lifecycle.

Odorant Receptors Expressing and Antennal Lobes Architecture Are Linked to Caste Dimorphism in Asian Honeybee, Apis cerana (Hymenoptera: Apidae).

Insects heavily rely on the olfactory system for food, mating, and predator evasion. However, the caste-related olfactory differences in Apis cerana, a eusocial insect, remain unclear. To explore the peripheral and primary center of the olfactory system link to the caste dimorphism in A. cerana, transcriptome and immunohistochemistry studies on the odorant receptors (ORs) and architecture of antennal lobes (ALs) were performed on different castes. Through transcriptomesis, we found more olfactory receptor genes in queens and workers than in drones, which were further validated by RT-qPCR, indicating caste dimorphism. Meanwhile, ALs structure, including volume, surface area, and the number of glomeruli, demonstrated a close association with caste dimorphism. Particularly, drones had more macroglomeruli possibly for pheromone recognition. Interestingly, we found that the number of ORs and glomeruli ratio was nearly 1:1. Also, the ORs expression distribution pattern was very similar to the distribution of glomeruli volume. Our results suggest the existence of concurrent plasticity in both the peripheral olfactory system and ALs among different castes of A. cerana, highlighting the role of the olfactory system in labor division in insects.

Range-wide genetic analysis of an endangered bumble bee (Bombus affinis, Hymenoptera: Apidae) reveals population structure, isolation by distance, and low colony abundance.

Declines in bumble bee species range and abundances are documented across multiple continents and have prompted the need for research to aid species recovery and conservation. The rusty patched bumble bee (Bombus affinis) is the first federally listed bumble bee species in North America. We conducted a range-wide population genetics study of B. affinis from across all extant conservation units to inform conservation efforts. To understand the species' vulnerability and help establish recovery targets, we examined population structure, patterns of genetic diversity, and population differentiation. Additionally, we conducted a site-level analysis of colony abundance to inform prioritizing areas for conservation, translocation, and other recovery actions. We find substantial evidence of population structuring along an east-to-west gradient. Putative populations show evidence of isolation by distance, high inbreeding coefficients, and a range-wide male diploidy rate of ~15%. Our results suggest the Appalachians represent a genetically distinct cluster with high levels of private alleles and substantial differentiation from the rest of the extant range. Site-level analyses suggest low colony abundance estimates for B. affinis compared to similar datasets of stable, co-occurring species. These results lend genetic support to trends from observational studies, suggesting that B. affinis has undergone a recent decline and exhibit substantial spatial structure. The low colony abundances observed here suggest caution in overinterpreting the stability of populations even where B. affinis is reliably detected interannually. These results help delineate informed management units, provide context for the potential risks of translocation programs, and help set clear recovery targets for this and other threatened bumble bee species.

Cloning, functional expression, and pharmacological characterization of inwardly rectifying potassium channels (Kir) from Apis mellifera.

Potassium channels belong to the super family of ion channels and play a fundamental role in cell excitability. Kir channels are potassium channels with an inwardly rectifying property. They play a role in setting the resting membrane potential of many excitable cells including neurons. Although putative Kir channel family genes can be found in the Apis mellifera genome, their functional expression, biophysical properties, and sensitivity to small molecules with insecticidal activity remain to be investigated. We cloned six Kir channel isoforms from Apis mellifera that derive from two Kir genes, AmKir1 and AmKir2, which are present in the Apis mellifera genome. We studied the tissue distribution, the electrophysiological and pharmacological characteristics of three isoforms that expressed functional currents (AmKir1.1, AmKir2.2, and AmKir2.3). AmKir1.1, AmKir2.2, and AmKir2.3 isoforms exhibited distinct characteristics when expressed in Xenopus oocytes. AmKir1.1 exhibited the largest potassium currents and was impermeable to cesium whereas AmKir2.2 and AmKir2.3 exhibited smaller currents but allowed cesium to permeate. AmKir1 exhibited faster opening kinetics than AmKir2. Pharmacological experiments revealed that both AmKir1.1 and AmKir2.2 are blocked by the divalent ion barium, with IC50 values of 10-5 and 10-6 M, respectively. The concentrations of VU041, a small molecule with insecticidal properties required to achieve a 50% current blockade for all three channels were higher than those needed to block Kir channels in other arthropods, such as the aphid Aphis gossypii and the mosquito Aedes aegypti. From this, we conclude that Apis mellifera AmKir channels exhibit lower sensitivity to VU041.

Population-wide modelling reveals prospects of marker-assisted selection for parasitic mite resistance in honey bees.

In 2019, a joint eight-variant model was published in which eight single nucleotide polymorphisms (SNPs) in seven Apis mellifera genes were associated with Varroa destructor drone brood resistance (DBR, i.e. mite non-reproduction in drone brood). As this model was derived from only one Darwinian Black Bee Box colony, it could not directly be applied on a population-overarching scale in the northern part of Belgium (Flanders), where beekeepers prefer the carnica subspecies. To determine whether these eight SNPs remained associated with the DBR trait on a Flemish colony-broad scope, we performed population-wide modelling through sampling of various A. mellifera carnica colonies, DBR scoring of Varroa-infested drone brood and variant genotyping. Novel eight-variant modelling was performed and the classification performance of the eight SNPs was evaluated. Besides, we built a reduced three-variant model retaining only three genetic variants and found that this model classified 76% of the phenotyped drones correctly. To examine the spread of beneficial alleles and predict the DBR probability distribution in Flanders, we determined the allelic frequencies of the three variants in 292 A. mellifera carnica queens. As such, this research reveals prospects of marker-assisted selection for Varroa drone brood resistance in honeybees.

Roles of DNA Methylation in Color Alternation of Eastern Honey Bees (Apis cerana) Induced by the Royal Jelly of Western Honey Bees (Apis mellifera).

Honey bees have a very interesting phenomenon where the larval diets of two different honey bee species are exchanged, resulting in altered phenotypes, namely, a honey bee nutritional crossbreed. This is a classical epigenetic process, but its underlying mechanisms remain unclear. This study aims to investigate the contribution of DNA methylation to the phenotypic alternation of a Apis mellifera-Apis cerana nutritional crossbreed. We used a full nutritional crossbreed technique to rear A. cerana queens by feeding their larvae with A. mellifera royal-jelly-based diets in an incubator. Subsequently, we compared genome-wide methylation sequencing, body color, GC ratio, and the DMRs between the nutritional crossbreed, A. cerana queens (NQs), and control, A. cerana queens (CQs). Our results showed that the NQ's body color shifted to yellow compared to the black control queens. Genome methylation sequencing revealed that NQs had a much higher ratio of mCG than that of CQs. A total of 1020 DMGs were identified, of which 20 DMGs were enriched into key pathways for melanin synthesis, including tryptophan, tyrosine, dopamine, and phenylalanine KEGG pathways. Three key differentially methylated genes [OGDH, ALDH(NAD+) and ALDH7] showed a clear, altered DNA methylation in multiple CpG islands in NQs compared to CQs. Consequently, these findings revealed that DNA methylation participates in A. cerana-A. mellifera nutritional crossbreeding as an important epigenetic modification. This study serves as a model of cross-kingdom epigenetic mechanisms in insect body color induced by environmental factors.

Phylogenomic analysis of Stylops reveals the evolutionary history of a Holarctic Strepsiptera radiation parasitizing wild bees.

Holarctic Stylops is the largest genus of the enigmatic insect order Strepsiptera, twisted winged parasites. Members of Stylops are obligate endoparasites of Andrena mining bees and exhibit extreme sexual dimorphism typical of Strepsiptera. So far, molecular studies on Stylops have focused on questions on species delimitation. Here, we utilize the power of whole genome sequencing to infer the phylogeny of this morphologically challenging genus from thousands of loci. We use a species tree method, concatenated maximum likelihood analysis and Bayesian analysis with a relaxed clock model to reconstruct the phylogeny of 46 Stylops species, estimate divergence times, evaluate topological consistency across methods and infer the root position. Furthermore, the biogeographical history and coevolutionary patterns with host species are assessed. All methods recovered a well resolved topology with close to all nodes maximally supported and only a handful of minor topological variations. Based on the result, we find that included species can be divided into 12 species groups, seven of them including only Palaearctic species, three Nearctic and two were geographically mixed. We find a strongly supported root position between a clade formed by the spreta, thwaitesi and gwynanae species groups and the remaining species and that the sister group of Stylops is Eurystylops or Eurystylops + Kinzelbachus. Our results indicate that Stylops originated in the Western Palaearctic or Western Palaearctic and Nearctic in the early Neogene or late Paleogene, with four independent dispersal events to the Nearctic. Cophylogenetic analyses indicate that the diversification of Stylops has been shaped by both significant coevolution with the mining bee hosts and host-shifting. The well resolved and strongly supported phylogeny will provide a valuable phylogenetic basis for further studies into the fascinating world of Strepsipterans.

Protein subcellular relocalization and function of duplicated flagellar calcium binding protein genes in honey bee trypanosomatid parasite.

The honey bee trypanosomatid parasite, Lotmaria passim, contains two genes that encode the flagellar calcium binding protein (FCaBP) through tandem duplication in its genome. FCaBPs localize in the flagellum and entire body membrane of L. passim through specific N-terminal sorting sequences. This finding suggests that this is an example of protein subcellular relocalization resulting from gene duplication, altering the intracellular localization of FCaBP. However, this phenomenon may not have occurred in Leishmania, as one or both of the duplicated genes have become pseudogenes. Multiple copies of the FCaBP gene are present in several Trypanosoma species and Leptomonas pyrrhocoris, indicating rapid evolution of this gene in trypanosomatid parasites. The N-terminal flagellar sorting sequence of L. passim FCaBP1 is in close proximity to the BBSome complex, while that of Trypanosoma brucei FCaBP does not direct GFP to the flagellum in L. passim. Deletion of the two FCaBP genes in L. passim affected growth and impaired flagellar morphogenesis and motility, but it did not impact host infection. Therefore, FCaBP represents a duplicated gene with a rapid evolutionary history that is essential for flagellar structure and function in a trypanosomatid parasite.

Physiological specialization of the brain in bumble bee castes: Roles of dopamine in mating-related behaviors in female bumble bees.

We aimed to investigate the roles of dopamine in regulating caste-specific behaviors in bumble bees and mating-related behaviors in bumble bee gynes. We examined caste differences in behaviors, biogenic amine levels, and expression levels of genes encoding dopamine receptors in the brains of bumble bees, and analyzed the effects of dopamine-related drugs on bumble bee behavior. Locomotor and flight activities were significantly higher in 8-day-old gynes and light avoidance was significantly lower in 4-8-day-old gynes than in same-aged workers. Brain levels of dopamine and octopamine were significantly higher in 8-day-old gynes than in same-aged workers, but tyramine and serotonin levels did not differ between the castes. Relative expression levels of the dopamine receptor gene BigDop1 were significantly lower in 8-day-old gynes than in same-aged workers, but expression levels of other dopamine receptor genes did not differ between castes. Dopamine significantly enhanced locomotor and flight activities in 7-9-day-old workers, whereas the dopamine receptor antagonist flupentixol inhibited flight activity and mating acceptance in same-aged gynes. These results suggest that dopamine plays important roles in gyne-specific behavior in bumble bees and has a common dopaminergic function in female eusocial bees.

Dynamics and regulatory role of circRNAs in Asian honey bee larvae following fungal infection.

Non-coding RNA (ncRNA) plays a vital part in the regulation of immune responses, growth, and development in plants and animals. Here, the identification, characteristic analysis, and molecular verification of circRNAs in Apis cerana cerana worker larval guts were conducted, followed by in-depth investigation of the expression pattern of larval circRNAs during Ascosphaera apis infection and exploration of the potential regulatory part of differentially expressed circRNAs (DEcircRNAs) in host immune responses. A total of 3178 circRNAs in the larval guts of A. c. cerana were identified, with a length distribution ranging from 15 to 96,007 nt. Additionally, 155, 95, and 86 DEcircRNAs were identified in the in the 4-, 5-, and 6-day-old larval guts following A. apis infection. These DEcircRNAs were predicted to target 29, 25, and 18 parental genes relevant to 12, 20, and 17 GO terms as well as 144, 114, and 61 KEGG pathways, including 5 cellular and 4 humoral immune pathways. Complex competing endogenous RNA (ceRNA) regulatory networks were detected as being formed among DEcircRNAs, DEmiRNAs, and DEmRNAs. The target DEmRNAs were engaged in 36, 47, and 47 GO terms as well as 331, 332, and 331 pathways, including 6 cellular and 6 humoral immune pathways. Further, 19 DEcircRNAs, 5 DEmiRNAs, and 3 mRNAs were included in the sub-networks relative to 3 antioxidant enzymes. Finally, back-splicing sites within 15 circRNAs and the difference in the 15 DEcircRNAs' expression between uninoculated and A. apis-inoculated larval guts were confirmed based on molecular methods. These findings not only enrich our understanding of bee host-fungal pathogen interactions but also lay a foundation for illuminating the mechanism underlying the DEcircRNA-mediated immune defense of A. c. cerana larvae against A. apis invasion. KEY POINTS: • The expression pattern of circRNAs was altered in the A. cerana worker larval guts following A. apis infection. • Back-splicing sites within 15 A. cerana circRNAs were verified using molecular approaches. DEcircRNAs potentially modulated immune responses and antioxidant enzymes in A. apis-challenged host guts.

Genomic characterization, phylogenetic and expression analysis of foraging gene in Apis mellifera.

The genomic characterization of the foraging gene and its expression analysis are required to better understand the behavior of honey bees (Apis mellifera). The present study performed a genome-wide characterization of the foraging gene, analyzing its physicochemical properties, phylogenetic features, and expression. An in silico analysis was carried out to characterize the foraging gene and the motifs and conserved domains of the encoded protein to predict its physicochemical properties. Moreover, a phylogenetic analysis of the foraging gene was performed in different species using MEGAX. The relative expression of the foraging gene was determined using qRT-PCR in two groups of forager bee samples (incoming and outgoing bees) during two seasons (five times per day). In addition, the queen effect was evaluated in another experiment. The results revealed that foraging gene expression and bee traffic were influenced by the interaction of season and daytime. The daily foraging traffic and transcription level of the foraging gene were the same in both seasons. The traffic of bees and the transcription abundance of the foraging gene were the highest in the middle and at the end of the day in the first and second seasons, respectively. Furthermore, the mRNA expression of the foraging gene was relatively higher in incoming bees than in outgoing bees. The queen also had a significant effect on the outgoing bees. We conclude that gene-environment interactions affect the foraging behavior of bees through the modulation of the foraging gene transcription.

Plant-pollinator network architecture does not impact intraspecific microbiome variability.

Variation in how individuals interact with food resources can directly impact, and be affected by, their microbial interactions due to the potential for transmission. The degree to which this transmission occurs, however, may depend on the structure of forager networks, which determine the community-scale transmission opportunities. In particular, how the community-scale opportunity for transfer balances individual-scale barriers to transmission is unclear. Examining the bee-flower and bee-microbial interactions of over 1000 individual bees, we tested (1) the degree to which individual floral visits predicted microbiome composition and (2) whether plant-bee networks with increased opportunity for microbial transmission homogenized the microbiomes of bees within that network. The pollen community composition carried by bees was associated with microbiome composition at some sites, suggesting that microbial transmission at flowers occurred. Contrary to our predictions, however, microbiome variability did not differ based on transfer opportunity: bee microbiomes in asymmetric networks with high opportunity for microbial transfer were similarly variable compared to microbiomes in networks with more evenly distributed links. These findings suggest that microbial transmission at flowers is frequent enough to be observed at the community level, but that community network structure did not substantially change the dynamics of this transmission, perhaps due to filtering processes in host guts.

Revealing the genome of the microsporidian Vairimorpha bombi, a potential driver of bumble bee declines in North America.

Pollinators are vital for food security and the maintenance of terrestrial ecosystems. Bumblebees are important pollinators across northern temperate, arctic, and alpine ecosystems, yet are in decline across the globe. Vairimorpha bombi is a parasite belonging to the fungal class Microsporidia that has been implicated in the rapid decline of bumblebees in North America, where it may be an emerging infectious disease. To investigate the evolutionary basis of pathogenicity of V. bombi, we sequenced and assembled its genome using Oxford Nanopore and Illumina technologies and performed phylogenetic and genomic evolutionary analyses. The genome assembly for V. bombi is 4.73 Mb, from which we predicted 1,870 protein-coding genes and 179 tRNA genes. The genome assembly has low repetitive content and low GC content. V. bombi's genome assembly is the smallest of the Vairimorpha and closely related Nosema genera, but larger than those found in the Encephalitozoon and Ordospora sister clades. Orthology and phylogenetic analysis revealed 18 core conserved single-copy microsporidian genes including the histone acetyltransferase (HAT) GCN5. Surprisingly, V. bombi was unique to the microsporidia in not encoding the second predicted HAT ESA1. The V. bombi genome assembly annotation included 265 unique genes (i.e. not predicted in other microsporidia genome assemblies), 20% of which encode a secretion signal, which is a significant enrichment. Intriguingly, of the 36 microsporidian genomes we analyzed, 26 also had a significant enrichment of secreted signals encoded by unique genes, ranging from 6 to 71% of those predicted genes. These results suggest that microsporidia are under selection to generate and purge diverse and unique genes encoding secreted proteins, potentially contributing to or facilitating infection of their diverse hosts. Furthermore, V. bombi has 5/7 conserved spore wall proteins (SWPs) with its closest relative V. ceranae (that primarily infects honeybees), while also uniquely encoding four additional SWPs. This gene class is thought to be essential for infection, providing both environmental protection and recognition and uptake into the host cell. Together, our results show that SWPs and unique genes encoding a secretion signal are rapidly evolving in the microsporidia, suggesting that they underpin key pathobiological traits including host specificity and pathogenicity.

How reliable is metabarcoding for pollen identification? An evaluation of different taxonomic assignment strategies by cross-validation.

Metabarcoding is a powerful tool, increasingly used in many disciplines of environmental sciences. However, to assign a taxon to a DNA sequence, bioinformaticians need to choose between different strategies or parameter values and these choices sometimes seem rather arbitrary. In this work, we present a case study on ITS2 and rbcL databases used to identify pollen collected by bees in Belgium. We blasted a random sample of sequences from the reference database against the remainder of the database using different strategies and compared the known taxonomy with the predicted one. This in silico cross-validation (CV) approach proved to be an easy yet powerful way to (1) assess the relative accuracy of taxonomic predictions, (2) define rules to discard dubious taxonomic assignments and (3) provide a more objective basis to choose the best strategy. We obtained the best results with the best blast hit (best bit score) rather than by selecting the majority taxon from the top 10 hits. The predictions were further improved by favouring the most frequent taxon among those with tied best bit scores. We obtained better results with databases containing the full sequences available on NCBI rather than restricting the sequences to the region amplified by the primers chosen in our study. Leaked CV showed that when the true sequence is present in the database, blast might still struggle to match the right taxon at the species level, particularly with rbcL. Classical 10-fold CV-where the true sequence is removed from the database-offers a different yet more realistic view of the true error rates. Taxonomic predictions with this approach worked well up to the genus level, particularly for ITS2 (5-7% of errors). Using a database containing only the local flora of Belgium did not improve the predictions up to the genus level for local species and made them worse for foreign species. At the species level, using a database containing exclusively local species improved the predictions for local species by ∼12% but the error rate remained rather high: 25% for ITS2 and 42% for rbcL. Foreign species performed worse even when using a world database (59-79% of errors). We used classification trees and GLMs to model the % of errors vs. identity and consensus scores and determine appropriate thresholds below which the taxonomic assignment should be discarded. This resulted in a significant reduction in prediction errors, but at the cost of a much higher proportion of unassigned sequences. Despite this stringent filtering, at least 1/5 sequences deemed suitable for species-level identification ultimately proved to be misidentified. An examination of the variability in prediction accuracy between plant families showed that rbcL outperformed ITS2 for only two of the 27 families examined, and that the % correct species-level assignments were much better for some families (e.g. 95% for Sapindaceae) than for others (e.g. 35% for Salicaceae).

Insights from Melipona bicolor hybrid genome assembly: a stingless bee genome with chromosome-level scaffold.

BACKGROUND: The highly eusocial stingless bees are crucial pollinators of native and agricultural ecosystems. Nevertheless, genomic studies within this bee tribe remain scarce. We present the genome assembly of the stingless bee Melipona bicolor. This bee is a remarkable exception to the typical single-queen colony structure, since in this species, multiple queens may coexist and share reproductive duties, resulting in genetically diverse colonies with weak kinship connections. As the only known genuinely polygynous bee, M. bicolor's genome provides a valuable resource for investigating sociality beyond kin selection. RESULTS: The genome was assembled employing a hybrid approach combining short and long reads, resulting in 241 contigs spanning 259 Mb (N50 of 6.2 Mb and 97.5% complete BUSCOs). Comparative analyses shed light on some evolutionary aspects of stingless bee genomics, including multiple chromosomal rearrangements in Melipona. Additionally, we explored the evolution of venom genes in M. bicolor and other stingless bees, revealing that, apart from two genes, the conserved repertoire of venom components remains under purifying selection in this clade. CONCLUSION: This study advances our understanding of stingless bee genomics, contributing to the conservation efforts of these vital pollinators and offering insights into the evolutionary mechanisms driving their unique adaptations.

The genome of the blind bee louse fly reveals deep convergences with its social host and illuminates Drosophila origins.

Social insects' nests harbor intruders known as inquilines,1 which are usually related to their hosts.2,3 However, distant non-social inquilines may also show convergences with their hosts,4,5 although the underlying genomic changes remain unclear. We analyzed the genome of the wingless and blind bee louse fly Braula coeca, an inquiline kleptoparasite of the western honey bee, Apis mellifera.6,7 Using large phylogenomic data, we confirmed recent accounts that the bee louse fly is a drosophilid8,9 and showed that it had likely evolved from a sap-breeder ancestor associated with honeydew and scale insects' wax. Unlike many parasites, the bee louse fly genome did not show significant erosion or strict reliance on an endosymbiont, likely due to a relatively recent age of inquilinism. However, we observed a horizontal transfer of a transposon and a striking parallel evolution in a set of gene families between the honey bee and the bee louse fly. Convergences included genes potentially involved in metabolism and immunity and the loss of nearly all bitter-tasting gustatory receptors, in agreement with life in a protective nest and a diet of honey, pollen, and beeswax. Vision and odorant receptor genes also exhibited rapid losses. Only genes whose orthologs in the closely related Drosophila melanogaster respond to honey bee pheromone components or floral aroma were retained, whereas the losses included orthologous receptors responsive to the anti-ovarian honey bee queen pheromones. Hence, deep genomic convergences can underlie major phenotypic transitions during the evolution of inquilinism between non-social parasites and their social hosts.

The adaptation of bumblebees to extremely high elevation associated with their gut microbiota.

Bumblebees are among the most abundant and important pollinators for sub-alpine and alpine flowering plant species in the Northern Hemisphere, but little is known about their adaptations to high elevations. In this article, we focused on two bumblebee species, Bombus friseanus and Bombus prshewalskyi, and their respective gut microbiota. The two species, distributed through the Hengduan Mountains of southwestern China, show species replacement at different elevations. We performed genome sequencing based on 20 worker bee samples of each species. Applying evolutionary population genetics and metagenomic approaches, we detected genes under selection and analyzed functional pathways between bumblebees and their gut microbes. We found clear genetic differentiation between the two host species and significant differences in their microbiota. Species replacement occurred in both hosts and their bacteria (Snodgrassella) with an increase in elevation. These extremely high-elevation bumblebees show evidence of positive selection related to diverse biological processes. Positively selected genes involved in host immune systems probably contributed to gut microbiota changes, while the butyrate generated by gut microbiota may influence both host energy metabolism and immune systems. This suggests a close association between the genomes of the host species and their microbiomes based on some degree of natural selection.IMPORTANCETwo closely related and dominant bumblebee species, distributed at different elevations through the Hengduan Mountains of southwestern China, showed a clear genomic signature of adaptation to elevation at the molecular level and significant differences in their respective microbiota. Species replacement occurred in both hosts and their bacteria (Snodgrassella) with an increase in elevation. Bumblebees' adaptations to higher elevations are closely associated with their gut microbiota through two biological processes: energy metabolism and immune response. Information allowing us to understand the adaptive mechanisms of species to extreme conditions is implicit if we are to conserve them as their environments change.